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Abstract 

It is shown, within the framework of the Ginzburg-Landau theory for a 

superconductor with d x 2_ y 2 symmetry, that the passing of a supercurrent 

through the sample results, in general, in the induction of order-parameter 

components of distinct symmetry. The induction of s-wave and d xy ( x 2_ y 2y 

wave components are considered in detail. It is shown that in both cases 

the order parameter remains gapless; however, the structure of the lines of 

nodes and the lobes of the order parameter are modified in distinct ways, and 

the magnitudes of these modifications differ in their dependence on the (a-b 

plane) current direction. The magnitude of the induced s-wave component is 

estimated using the results of the calculations of Ren et al. [Phys. Rev. Lett. 

74, 3680 (1995)], which are based on a microscopic approach. 
PACS numbers: 74.20.De,74.25.Fy,74.72.-h,74.76.Bz 
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I. INTRODUCTION 



As a result of recent experimental and theoretical [ffl work there is an emerging 

consensus that the symmetry of the superconducting state in the high-temperature super- 
conducting materials is that of d x 2_ y 2-wave. Given this situation, it seems worthwhile to 
explore the phenomenological implications of such a state, even if the microscopic origin of 
the superconductivity has not yet been fully established. In particular, the following fea- 
ture of the phenomenological theory of d-wave superconductors has attracted the attention 
of a number of workers. In the absence of any external agents (e.g., magnetic fields, sur- 
faces, currents, etc.), the only component of the superconducting order parameter that has 
a non-zero mean equilibrium value is the component with d x 2_ y 2 symmetry, the other (i.e. 
subsidiary) components exhibiting equilibrium fluctuations around a mean value of zero. 
External forces can give rise to non-zero mean values of subsidiary components. Interest- 
ingly, general symmetry considerations permit a coupling between the gradients of d x 2_ y 2 
and of other components ||. In particular, this means that any inhomogeneity in d x 2_ y 2 
acts as a source of inhomogeneities in other components and, therefore, as a source of these 
components themselves. This mechanism has been exploited by a number of authors. No- 
tably, Volovik ||, followed by Soininen et al. and other workers p3| , |TT|| , have predicted 
that the vortices in a c?-wave superconductor should exhibit a rich structure, in which s- and 
o?-components of the order parameter co-exist. Furthermore, the surface regions of these 



superconductors are predicted to be in a mixed s-rf-state [|13 



In this work we pursue one further consequence of the gradient coupling mechanism, 
viz., we predict that an external current can induce non-zero subsidiary components of the 
superconducting order parameter via the (current-induced) inhomogeneity of the dominant 
(i.e. d x 2_ y 2) component. As opposed to the cases of vortices and surfaces |L3 



which both have an amplitude variation of the dominant {d x 2_ y 2) component, the induction 
of subsidiary components by the current requires only its phase variation. In Section II, 
we present the Ginzburg-Landau theory of the current-induced s-component in a <i-wave 
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superconductor. Our treatment is based on the Ginzburg-Landau (GL) theory for a <i-wave 
superconductor [P,[n]l that incorporates the effects of s/ci-coupling. As it is not yet clear 
which of the subsidiary components has the strongest coupling to the d x 2_ y 2 -component, we 
then (in Section III) extend this treatment to include subsidiary irreducible representations 
of the tetragonal (-D4) symmetry group other than s-wave. Finally, we discuss some ex- 
perimental settings in which current-induced subsidiary components of the order parameter 
might be observable. 



II. THE CASE OF s/d x 2_ y 2 COUPLING 

First we focus on the case of sjd x i_ y i -coupling, in which the order parameter has two 
spatially varying complex components s(r) and d(r). We neglect the magnetic fields induced 
by the current jl2|| . The GL equations for d- and s-components were derived in Refs. [f>l,|TTl; 
they are 

(-7rfV 2 + a d )d + 7 ,(V 2 - V 2 )s + 2(3 2 \d\ 2 d + (3 3 \s\ 2 d + 2[3 4 s 2 d* = 0, (la) 
(-7 S V 2 + a s )s + lv (V 2 x - V 2 y )d + 2(3 1 \s\ 2 s + (3 3 \d\ 2 s + 2{3 4 d 2 s* = 0, (lb) 

where 7 P = h 2 /2m p , and p — d,s, v. The current density is given by 

J = —{d*Vd - c.c.} + — {s*Vs - c.c.} 
irrid im s 

- x {s*V x d - dV x s* - c.c.} + y {s*V y d - dV y s* - c.c.}, (2) 

im v im v 

where we have chosen the effective charge e to be twice the electron charge. The parameters 
of the GL equations (jT^,[T5|) are chosen in such a way ||[lT| that, in the absence of the 



current, \d\ > and s = 0. In the presence of the current, we assume that s is nonzero 
but small compared with d and, therefore, can be analyzed perturbatively. In this way, the 
inhomogeneity in d acts as a source for s. To the zeroth order in perturbation theory, we 
have 

d = d e i(lor , s = s = 0. (3) 



The amplitude do and the wave vector qo are found in the usual way |14| from Eqs. fllajpl) 
with s having been set to zero: 

Qo = C *(1 _ f 2 ) 1/2 , and qo || J labelEQ : q, (4a) 

j = J/J c = 3VSf 2 (l-f) 1/2 /2, (4b) 



where / = do/J\ctd 1/2/32 is the dimensionless (i-wave order parameter normalized by its 



equilibrium value, = J fi 2 /2md\ad\ is the correlation length of the <i-wave order parameter, 
and J c is the critical current density. The dependence of /, and thus of do, on j is given 
by the implicit relation (|4b]) . In particular, / = 1 for j = 0, and / approaches the value of 



2/3 from above as j approaches 1 from below; for j > 1 we have / = 0. To first order in 
perturbation theory, s acquires a non-zero value and d changes from its zeroth-order value. 
This also leads to a change in the right-hand side of Eq. (pj), which determines the wave 
vector of the order parameter for a given current density. This means that at this order the 
wave vector found in the previous order is changed. Thus the appropriate Ansatz at first 
order is 

d = d e i(qo+qi) - r + d 1 e < ( qo+qi > ,r , (5a) 
8 = Sl e i(qo+qi) - r , (5b) 

where quantities with the subscript 1 are small compared to those with the subscript 0. As 
can be readily checked, this Ansatz satisfies Eqs. ([Tl^,|iT)|,|2|) . Keeping only terms linear in di, 
S\ and qi, and after some lengthy but straightforward algebra ||15|| , we obtain: 
£i = IvQo cos(20) 

* 7.rf + a. - ZS-Zl l SU - ^ + (ft + w' 

where is the angle between J and the x-axis in the a-6-plane. Already at this stage 
two conclusions can be made. First, the amplitude of the induced s-component depends 
not only on the amplitude but also on the direction of the current: \s\ is maximal for a 
current flowing along the major crystallographic axes (i.e., for = Oor0 = 7r/2) and is 
zero (at this order of perturbation theory) for a current flowing along the diagonal of the 



unit cell (i.e., for <fi = 7r/4). Second, the rather cumbersome expression @ is simplified 
considerably for temperatures T very close to the critical temperature T c (i.e., in the critical 
regime, when the GL approach is strictly valid). In the limit T —>■ T c all the terms in the 
denominator of Eq. (||) that contain go become small because diverges, and the last term 
in the denominator is small because do is small. On the other hand, as T c is not a critical 
temperature for the s-wave component, a s is nonzero in this limit; thus a s dominates the 
denominator. Equation (P) then takes on the simpler form: 

Si m d \a d \ 2 

— = 1 - / cos 20 . 7 

d m v a s 

Equation (]?]) shows that for generic values of the parameters, i.e., for ~ m v and 0—1, 
the smallness of Si with respect to d , and thus the validity of the perturbation theory, is 
guaranteed by the smallness of the ratio |a^|/a s oc (T c — T). Therefore, in the critical region, 
the perturbation theory is valid even for currents that are not small compared to J c (i.e., 
for values of / that are not close to 1). 

In order to obtain (semi-)quantitative estimates for the amplitude of the induced s-wave 
order parameter as given by Eq. (Q), we need to know the values of the phenomenological 
parameters of the GL theory. These can be estimated, e.g., by comparing Eqs. (la, |TFJ^) with 



the GL equations that were derived recently by Ren et al. [JTQ] from the Gor'kov equations for 
a particular microscopic model of pairing interactions fl6j . Reference flO| gives the following 
ratios of the phenomenological GL parameters: 

m s : md '■ m v = 1:2:2, (8a) 
Pi : fa ■ P 3 : Pa = 1 : 3/8 : 2 : 1/2, (8b) 
\a d \/a s = A d ln(T c /T)/2(l + V s /V d ), (8c) 

where is the BCS coupling constant in the ci-wave channel, and V d and V s are interaction 
parameters, which in the model of Ref. [0] describe nearest-neighbor attraction and on- 



site repulsion, respectively. As mentioned in Ref. |1U||, the s-wave component is induced by 



inhomogeneities in the ci-wave component even if V s = 0. For lack of better information 



about V s , we set it to zero, which does not significantly affect our results. By using the 
ratios of the GL parameters given above, Eq. (||) is cast into the following form: 

si 1 (1 - f 2 ) cos(20) 



d 2 1 9 (1-/*)W(20) 

A d ln(T c /T) 4 3/2-1 + i 



(9) 



To estimate the BCS coupling constant A^, we use the result of Monthoux and Pines [|T7 
who find that, in a spin-fluctuation model, the value of T c = 90 K is obtained for A^ close 
to 1 (the precise value depending on the doping). Solely for illustrative purposes, we use 
the value A^ = 1. The dependence of si/d® on j is shown in Fig. [I] for three values of 
t = (T c — T)/T c . For the values t = 0.01 and 0.1 the result given by Eq. (|j) is very close 
to that given by its simplified version Eq. (|7|). For t = 0.5 the result given by Eq. (|9|) is 
approximately one half of that given by Eq. (|7|). At temperatures far below T c (i.e. for 
t ~ 1), the microscopic derivation of the GL parameters leading to Eq. (U) is not strictly 
valid, and the term ln(T c /T) is expected to be replaced by t (note that ln(T c /T) = t for 
i<l) , thus avoiding the apparent singular behavior in Eq. (|9|). 

The presence of the s-wave component, which according to Eqs. ( |5l5||5"15| ) is in-phase 
with the d-wave component, implies that excitations with momentum along the $ = ±7r/4 
directions are no longer gapless. Rather, they have the energy gap given by A s = A d si/d , 
where is the maximum value of the <i-wave gap in the absence of the currents. The 
lines of nodes, oriented along the $ = ±7r/4 directions in the absence of current, are now 
rotated by the angle 5$ = \ cos _1 (A s /Ad) [see Figs. (@a,b)]. The fc-space structure of 
the order parameter undergoes an orthorhombic distortion, i.e., the current-induced s-wave 
component mimics the effect of having an orthorhombic (rather than tetragonal) lattice and 
no supercurrent. By using the typical value of A^ ~ 100 K, we see, e.g., that for t = 0.5 
(i.e., for T = 0.5T C ) and for j = 0.5 the gap A s ~ IK, and the lines of nodes rotate by 
5$ ~ 0.3°. We must emphasize that at lower temperatures (i.e., when t ~ 1) and for 
currents comparable to the critical current, there is no longer a natural small parameter in 
the theory that would automatically guarantee the smallness of Si with respect to d$. The 
fact that s\ remains small even in this region is due to the particular choice of the ratios of 

6 



the GL parameters. However, the GL theory is not expected to be quantitatively correct in 
this region, so the microscopic theory might give other numerical values of A s and 5$. The 
absence of a natural small parameter suggests that these values might be larger than those 
given by the perturbative treatment of the GL equations. 



III. COUPLING TO OTHER SUBSIDIARY ORDER-PARAMETER 

COMPONENTS 

So far, we have considered the coupling of the dominant d^.^-component to the (sub- 
sidiary) s-wave component, which is taken as the main subsidiary component in the mi- 
croscopic approaches of f^jllj. In general, the GL theory should incorporate all irreducible 
representations of the D± symmetry group; it is then the task of a microscopic theory to 
determine the dominant, and leading subdominant, components. Although the growing con- 
sensus is that the leading component corresponds to the d x 2_ y 2 representation, it is not clear, 
at present, which representation describes the subleading component ||13|| . Therefore, we now 



extend the treatment presented above to include couplings between the d x 2_ y 2 component 
and components of the order parameter other than s-wave. 

The irreducible representations of the (planar) D4 group are (see, e.g., Ref. 0): Ti or s- 
wave (transforming as x 2 +y 2 = 1), r% [transforming as xy(x 2 —y 2 )], or d x 2_ y 2 -wave (trans- 
forming as x 2 — y 2 ), (transforming as xy), and (transforming as the two-component 
vector {x, y}). (These representations are also commonly denoted A ig , A 2g , Bi g , B 2g , and 
E 2g , respectively.) Note that T 5 is a two-dimensional representation, whereas the other 
representations are one-dimensional. 

We now focus on the determination of the terms in the GL free energy describing the 
couplings between the gradients of T 3 {d x 2_ y 2) and of other representations. We consider 
only the leading terms of this type, i.e., terms of the form: 

CVV^V^r,, (10) 
where ipr k is the component of the order parameter transforming according to representation 



rV Here, /i, v = x, y, and i = 1, . . . , 5. These terms transform as the (reducible) representa- 
tion T = T 3 x Tj x T 5 x T 5 . As each term in the free energy must transform as a scalar, the 
maximum number of such gradient-coupling terms is given by the number Ni of times the 
identity (Ti) representation occurs in the decomposition of T into the irreducible represen- 
tations 0. Ni is given by the normalized product of characters corresponding to irreducible 
representations T 3 and Tk (see, e.g., Ref. This gives: Ni = 1 for % = 1, . . . , 4, and 

As = 0. First, we consider the case i — 2. A term satisfying all the symmetries of the group 
_D 4 can be written as 

C 

-{d x ipT 2 dyipr 3 + 9 y ipr 2 d x ipr 3 } + cc, (11) 

and, as N 2 = 1, there are no further independent terms 0. Next, we consider % — 4. The 
symmetry {x — > y,y — > x} imposes the conditions: C xx = —C yy and C xy = C yx = 0, 
while, e.g., the symmetry {x — > —x,y — > y} requires that C xx = 0. Therefore, all the 
constants are zero and there no are gradient-coupling terms to leading order for i — 4. The 
analysis of cases i — 1, 3 has been performed in Ref. f|, leading to Eqs. flla| , [TED , and N 5 is 



zero. Thus, the only case for which the induction of the subsidiary component of the order 
parameter by the current remains to be considered is that of the coupling between d x 2_ y 2 
and r 2 -representation [the latter is henceforth being referred to as d xy r x 2_ y z\\. 

We denote the component of the order parameter corresponding to the d X y^ x 2_ y 2^ - 
representation by a(r). In order to construct the GL free energy for the case of d xy ( x 2_ y 2y 
d x 2_ y 2 coupling, we: (i) note that the structure of terms other then the mixed gradient 
terms is the same as for the case of the s/d x 2_ y 2 coupling; and (ii) make use of Eq. (]TlD for 
the mixed gradient term. The usual variational procedure then leads to the following GL 
equations for d and a: 

(~7dV 2 + a d )d - -f w V x V y a + 2a 2 \d\ 2 d + a 3 \a\ 2 d + 2a 4 a 2 d* = 0, (12a) 
(~7a V 2 + a a )a - -f w V x V y d + 2ai|a| 2 a + a 3 \d\ 2 a + 2a 4 d 2 a* = 0, (12b) 

where 7 P = h 2 /2m p , and p = d,a,w. The current density takes the form 



-{d*Vd - c.c.} + - — {a*Va - c.c.} 



+ x- {a*V y d — c.c.} + y- {a*V x d — c.c.}. (13) 

[As the two last terms in Eq. ([13]) come from the variation of the (covariant) mixed gradient 

terms in the free energy with respect to the vector-potential, their structure is different from 

that of the analogous terms in Eq. @]. As in the case of s/d x 2_ y 2 -coupling, we assume 

that the amplitude of a induced by the current is small compared to d. The first-order 

perturbative calculation analogous to that for the case of the s/d x 2_ y 2 -coupling leads to the 

following result for the induced d xy ( x 2_ y 2- ) -component a,\\ 

<h = h^o sin(20) 

d ° ™$ +<*« + § Xi^2* + ffi° 2 + ( g 3 + 
We see that in contrast to the case of the s/d x 2_ y 2 -coupling [cf. Eq. (|J)], the induced 

dxy(x 2 -y 2 )~ component is zero for currents flowing along the principal crystallographic axes 

in the a — b plane (i.e., for = Oor0 = 7r/2) and reaches its maximum absolute value 

for currents flowing along the diagonal of the unit cell (i. e., for = vr/4). This difference 

could be used in an experiment to determine which of the two couplings (i.e., s/d x 2_ y 2 or 

d xy (x 2 -y 2 )/ 'd x 2_ y 2 ) is realized in a given HTS material. In the critical region, i.e., when 

T — > T c , the term a a dominates the denominator of Eq. fl!4D , due to the reasons described 

in the discussion of the s/d x 2_ y 2 -coupling, and Eq. ([14]) then takes the simpler form: 

— = -- (1 - / sm(20). (15) 

d 2m w a a 

As we are not aware of any microscopic theory describing the case of the d xy r x 2_ y 2) / d x 2_ y 2 
-coupling, we do not know the values of the GL parameters in Eqs. ( |12aj , |12b|) and, therefore, 
cannot give a quantitative estimate for the amplitude of the induced order parameter. 



IV. DISCUSSION AND CONCLUSIONS 



We have seen that the current-induced s-wave component introduces an orthorhombic- 
like distortion of the A;-space structure of the order parameter [Figs. (Qa,b)]. In contrast, the 
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induced d^^^-wave component distorts the fc-space structure as indicated in Fig. (0c). 
Note that the lines of nodes at $ = ±7r/4 do not rotate in the d xy ( x 2_ y 2ywave case, and, 
provided that a\ < 2do, no new nodes are introduced. The distortion of the zero-current, 
tetragonal structure of Fig. (fja) to the structure of Fig. (^|b) or Fig. (fjc) (or to a mixture 
of the latter two) by an externally imposed current may be experimentally observable using 
directional probes of the order parameter. Techniques such as photoemission or tunneling 
may be appropriate, provided that sufficiently high resolution can be obtained. 
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FIGURES 



0.05 




FIG. 1. The magnitude si of the s-wave component induced by the applied current density J, 
for the reduced temperature values t = 0.5 (dotted line), t = 0.1 (dashed line), and t = 0.01 (solid 
line). The curves are obtained from Eqs. (4b) and (9) with <\> = 0. 
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Fig. 2a 



Fig. 2b 



Fig. 2c 



FIG. 2. (a) The /e-space structure of the superconducting order parameter with: (a) pure 
d x 2_ y 2 symmetry [cos(2</>)] ; (b) mixed d x i_ y i and s symmetry [cos(2</>) + 0.1] ; (c) mixed d x 2_ y 2 
and d xy ( x 2_ y 2} symmetry [cos(2c/>) + O.3sin(20) cos(2</>)]. Note that for the purpose of illustration, 
the magnitude of the s-wave component in Fig. (b) has been chosen to be larger than that expected 
from microscopic estimates (see the main text). 
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